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Abstract 



O 

O ■ Nowadays there is no universally accepted definition of quantum chaos. 

-^ , In this paper we review and critically discuss different approaches to the 

XI 

Vh ' subject, such as Quantum Chaology and the Random Matrix Theory. Then 

we analyze the problem of dynamical chaos and the time scales associated 
with chaos suppression in quantum mechanics. 
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1 Introduction 

The aim of this paper is to review and discuss various definitions and ap- 
proaches to quantum chaos}^^ 

As yet, there is no universally accepted definition of quantum chaos. On 
the contrary, the meaning of classical chaos is beyond question. In classical 
mechanics a trajectory z(t) in the phase space Vt is chaotic if its maximal 
Lyapunov exponent A is positive. The Lyapunov exponent is defined as 

\=\im-\n\uj{t)\, (1) 

where uj{t) is a tangent vector to z(t) with the condition that |c(j(0)| = 1. The 
exponential instability of chaotic trajectories implies a continuous frequency 
spectrum of motion. The continuous spectrum, in turn, implies correlation 
decay; this property, which is called mixing in ergodic theory, is the most 
important property of dynamical motion for the validity of the statistical 
description.^"^ 

The problem of quantum chaos arose because the above mentioned con- 
dition of continuous spectrum for classical chaos is violated in quantum me- 
chanics. Indeed the energy and the frequency spectrum of any quantum mo- 
tion, bounded in phase space, are always discrete due to the non-commutative 
geometry (discreteness) of the phase space. According to the theory of dy- 
namical systems, such motion corresponds to the limiting case of regular 
motion. It means that there is no classical-like chaos at all in quantum 
mechanics. ^~^ 

Nevertheless, we shall show that it is reasonable and useful to apply 



the word chaos also in quantum mechanics. In the first part of this article 
some definitions of quantum chaos for stationary systems are discussed, while 
in the second part, the time evolution of classical and quantum systems is 
compared. 

2 Quantum Chaology and Spectral Statistics 

As is well known, in the study of the transition from order to chaos in classical 
systems a useful tool is the examination of the phase space properties, such 
as the Poincare sections.^ Such plots are not directly available in the case of 
quantum systems. In many papers, the Berry definition of quantum chaos is 
adopted: " Quantum Chaology is the study of semiclassical, but not classical, 
behaviour characteristic of systems whose classical motion exhibits chaos" .^ 

The idea, also suggested by authors like Percival and Gutzwiller,^ is to 
connect the behaviour of the eigenvalues and eigenfunctions of a quantum 
system to the different structure of the phase space of the corresponding 
classical system in the regular and chaotic region. 

In the context of quantum chaology, the spectral statistics of the energy 
levels are of great importance. Mehta defined: "A spectral statistic is a 
quantity which can be calculated from an observed sequence of levels alone, 
without other information and whose average value and variance are known 
from the theoretical model. A suitable statistic is one which is sensitive for 
the properties to be compared or distinguished and is insensitive for other 
details".^ In particular, it has been found that the spectral statistics of sys- 
tems with underlying classical chaotic behaviour and time-reversal symmetry 



agree with the predictions of the Gaussian Orthogonal Ensemble (GOE) of 
Random Matrix Theory (RMT),^ whereas quantum analogs of classically in- 
tegrable systems display the characteristics of the Poisson statistics. ^"^ Note 
that if the chaotic system is without time-reversal symmetry then, instead 
of the GOE, it follows the predictions of the Gaussian Unitary Ensemble 
(GUE). 

The most used spectral statistics of the energy levels are P{s) and A3(L). 
P{s) is the distribution of nearest-neighbour spacings Sj = (-Ej+i — Ei) of the 
unfolded levels Ei. It is obtained by accumulating the number of spacings 
that lie within the bin (s, s + As) and then normalizing P{s) to unit. As 
shown by Berry,^"'-*^^ for quantum systems whose classical analogs are inte- 
grable, P{s) is expected to follow the Poisson distribution 

P{s) = exp (-s) . (2) 

On the other hand, quantal analogs of chaotic systems exhibit the spectral 
properties of GOE with 

P(s) = ^sexp(-p), (3) 

which is the so-called Wigner distribution. Note that for systems without 
time-reversal symmetry the GUE predicts P{s) = (32/7r^)s^exp (— 4s^/7r). 

The statistic A3(L) is defined, for a fixed interval (— L/2, L/2), as the 
least-square deviation of the staircase function N{E) from the best straight 
line fitting it: 

AgfL) = - min / \N(E) - AE - B?dE , 

^ ' L A,B J^L/2- ^ 



where N{E) is the number of levels between E and zero for positive energy, 
between —E and zero for negative energy. The A3(L) statistic provides a 
measure of the degree of rigidity of the spectrum: for a given interval L, 
the smaller A3(L) is, the stronger is the rigidity, signifying the long-range 
correlations between levels. For this statistic the Poissonian prediction is 

ML) = ^ . (4) 

The GOE predicts the same behaviour for L << 1; instead for L >> 1 it 
gives 

A3(L) = ^logL. (5) 

In the GUE case one has A3(L) = (l/27r^) logL. It is useful to remember 
that Berry^^ has shown that A3(L) deviates from the universal predictions 
of RMT for very large L. 

Another probe, which is generally regarded (see for instance Ref. 14) as 
very sensitive to the structure of chaotic states, is the transition probability. 
For reasons of space, we mention only the pioneering work of French and 
his coworkers^^ and two more recent works concerning the interacting-boson 
model^^ (IBM) and the three-level Lipkin, Meshkov, Glick (LMG) model. ^^ 
In both works the results based on the transition probabilities between eigen- 
states of the system completely agree with the spectral statistics P{s) and 
A3(L). 

It is important to stress that even though the classical system is not 
known, to distinguish between ordered and chaotic states, the spectral statis- 
tics and the transition probabilities can be used.^^~^^ 



3 From Poisson to GOE Transition: Compar- 
ison with Experimental Data 

The agreement between the classical order-chaos transition and the quantal 
Poisson-GOE transition has been tested in many theoretical models, rang- 
ing from simple billiards^*^"^^ to more realistic systems like nuclei^^"^^ and 
elementary particles. ^^^^^ 

In this section we shall compare the transition Poisson-GOE with the ex- 
perimental data of two different systems: the atomic nuclei and the Hydrogen 
atom in a static magnetic field. 

3.1 Atomic Nuclei 

In atomic nuclei, as in other many-body systems, ordered and chaotic states 
generally coexist. '^^ 

a) The Low Energy Region 

The behaviour of spectral statistics near the ground state has been studied 
by Garret, German, Courtney and Espino^^ and Shriner, Mitchell and Von 
Egidy^^. The main results of these authors have been shown in Figure 1 
and Figure 2. As can be seen from the figures, in the low energy region the 
spectral statistics are in agreement with the Poisson ensemble or intermediate 
between Poisson and GOE. 

b) The High Energy Region 



Neutron resonance spectroscopy on a heavy even-even nucleus typically 
leads to the identification of about 150 to 170 s-wave resonances with J"^ = ^ 
located 8-10 MeV above the ground state of the compound system, with 
average spacings around 10 eV and average total widths around 1 eV. Proton 
resonance spectroscopy yields somewhat shorter sequences of levels with fixed 
spin and parity, with typically 60 to 80 members. 

For the statistical analysis, it is essential that the sequences be pure (no 
admixture of levels with different spin or parity) and complete (no missing 
levels). Only such sequences were considered by Haq, Pandey and Bohigas.^'^ 
Scaling each sequence to the same average level spacing and lumping together 
all sequences one leads to the "Nuclear Data Ensemble" (NDE), which con- 
tains 1726 level spacings. 

As shown in Figure 3, the agreement between the experimental data and 
the GOE predictions is surprisingly good (in the GOE model there are no 
free parameters). 

3.2 The Hydrogen Atom in the Strong Magnetic Field 

We now discuss the local statistical properties of energy levels of a Hydrogen 
atom in a uniform strong magnetic field. The Hamiltonian of the system is 
given by 

H = f-^--fL.^'l^i.'^y'), (6) 

2m r 2m Sm 

where the magnetic field B breaks the time-reversal symmetry. Although Eq. 
(6) is not time-reversal invariant, it can easily be written in a time-reversal 
invariant form.^^ In fact, the paramagnetic interaction ^Lz simply shifts 
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the whole series of levels with a fixed quantum number M (eignenstate of 
Lz with eigenvalue MK) and can be taken into account in the standard way: 
the Zeeman effect. 

The Hamiltonian (6) written in atomic units m = \q\ = 47reo = h = 1 is 

^ = y-^ + y(-^ + ^^), (7) 

where 7 = B /B^ is the magnetic field in atomic units and B^ = 2.35 ■ 10^ 
Tesla. This equation can be numerically solved for different values of the 
scaled energy e = E/(2'jY^^. Once the eigenvalues have been obtained, the 
spectral statistics P{s) and A3(L) can be calculated. Figure 4 shows the 
function P{s) for different values of the scaled energy e. Increasing e, a 
smooth Poisson-GOE transition can be observed. Figure 5 shows the spec- 
tral rigidity A3(L) in three different energy intervals. For this statistic the 
transition Poisson-GOE is also very clear. 

In addition, a comparison has been made between the theoretical energy 
levels and the experimental ones; the agreement is excellent. 

4 Quantum Chaos and Field Theory 

In the last few years there has been much interest in chaos in field theories. 
It is now well known that the spatially uniform limits of scalar electrodynam- 
ics and Yang-Mills theory exhibit classical chaotic motion. ^^ In this section 
we discuss quantum chaos in a field-theory schematic model, namely the 
spatially homogeneous SU(2) Yang-Mills-Higgs (YMH) system.^^^^^ The La- 



grangian density of the SU(2) YMH system is given by 



i(Z^,0)+(D»-n0)-i_ 



L = -{D,(l>y{Dy) - V{4>) - -F;^F^^'' , (8) 



where 



{D,<p) = d,<p-igA';r'<t>, (9) 

F;. = d,K - duA; + ge^'^A'^Al , (10) 

with T* = o''^/2, b = 1,2,3, generators of the SU(2) algebra, and where the 
potential of the scalar field (the Higgs field) is 

\/(0)=/i2|0p+Ai0r (11) 

In the (2+l)-dimensional Minkowski space {fi = 0,1,2) and with spatially 
homogeneous Yang-Mills and the Higgs fields 

d,A'; = d,<j) = 0, ^ = 1,2 (12) 

one considers the system in the region in which space fluctuations of fields 
are negligible compared to their time fluctuations. 

In the gauge Aq = and using the real triplet representation for the Higgs 
field one obtains 

L='^ + 1{AI + AI) - g'[\AlAl - \{A, . A,f+ 

+ {Al + Al)$' - {A, ■ $f - (I2 ■ 0)'] - V(0) , (13) 

where = {(f)\ (j)^ , (f)^) , Ai = {A\,AlA^) and I2 = (AlAlA^). 

When /i^ > 0, the potential V has a minimum at |0| =0, but for fi^ < 
the minimum is at 



which is the non zero Higgs vacuum. This vacuum is degenerate, and after 
spontaneous symmetry breaking the physical vacuum can be chosen 0o = 
(0,0,f). If A\ = qi, A2 = q2 and the other components of the Yang-Mills 
fields are zero, in the Higgs vacuum the Hamiltonian of the system reads 

H=l{pl+ pD + g'v'{ql + ql) + \g\lql , (14) 

where pi = gi and P2 = Q2- Here w'^ = 2g'^v'^ is the mass term of the Yang- 
Mills fields. This YMH Hamiltonian is a toy model for classical non-linear 
dynamics, with the attractive feature that the model emerges from particle 
physics. At low energy the motion near the minimum of the potential 

V{qi, gs) = 9'v\ql + ql) + \g\lql , (15) 

where the Gaussian curvature is positive, is periodic or quasiperiodic and is 
separated from the instability region by a line of zero curvature; if the energy 
is increased, the system will be for some initial conditions in a region of 
negative curvature, where the motion is chaotic. According to this scenario, 
the energy E^. of chaos-order transition is equal to the minimum value of the 
line of zero Gaussian curvature K{qi,q2) on the potential-energy surface. It 
is easy to show that the minimal energy on the zero-curvature line is given 
by: 

Ec = VrrUK = 0, gi) = 6g\^ , (16) 

and by inverting this equation one obtains Vc = {E/6g'^Y^*. There is an 
order-chaos transition by increasing the energy E of the system and a chaos- 
order transition by increasing the value v of the Higgs field in the vacuum. 
Thus, there is only one transition regulated by the sole parameter E/{g'^v'^). 
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It is important to point out that in general the curvature criterion guar- 
antees only a local instabilit'i/^^ and should therefore be combined with the 
Poincare sections^^-* . Chaotic regions on the surface of the section are char- 
acterized by a set of randomly distributed points, and regular regions by 
dotted or solid curves. Figure 6 shows the Poincare sections, which confirm 
the analytical predictions of the curvature criterion: the critical value of the 
onset of chaos is in very good agreement with the Poincare sections. 

In quantum mechanics the generalized coordinates of the YMH system 
satisfy the usual commutation rules [qk,Pi] = i^ki, with k,l = 1,2. Introduc- 
ing the creation and destruction operators 

the quantum YMH Hamiltonian can be written 

H = Ho + \g'y , (18) 

where 

Hq = uj{a'lai + 0,20,2 + 1) , (19) 

V = 2^{ai + at)\a2 + ot)\ (20) 

with cj^ = 2g'^v'^ and [a^, a^] = Ski, k,l = 1, 2. If \n1n2 > is the basis of the 
occupation numbers of the two harmonic oscillators, the matrix elements are 

< n[n2\Ho\nin2 >= uj{ni + ns + l)5„;„/„;„2 ' (^1) 

and 

< n[n'2\V\n,n2 >= —[^ni{ni - 1)'5„;„,_2+V(^i + l)(^i + '^K[m+2+i^ni+'^K[n^^ 
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X [v/n2(n2-l)<5„,„^_2 + ^K + l)(n2 + 2)5„.„^^, + (2n, + 1)5„.„J . (22) 

Figure 7 shows the P(s) distribution for different values of the parameter v. 
The figure shows a Wigner-Poisson transition by increasing the value v of 
the Higgs field in the vacuum. The P{s) distribution is fitted by the Brody 
function 

P(s,cu) =a(cj + l)s'^exp(-as'^+i) , (23) 

with 

a = {V[^]r' . (24) 

This function interpolates between the Poisson distribution {uo = 0) of inte- 
grable systems and the Wigner distribution {uj = 1) of chaotic ones, and thus 
the parameter u can be used as a simple quantitative measure of the degree 
of chaoticity. By using the P(s) distribution and the Brody function it is 
possible to give a quantitative measure of the degree of quantal chaoticity of 
the system. The numerical calculations of Figure 6 and 7 clearly show the 
quantum chaos-order transition and its connection to the classical one. 

5 Alternative Approaches to Quantum Chaos 

A different approach to quantum chaos has been discussed by Sakata and 
his coworkers. ^^ The example discussed by Sakata is a system of an even 
number of fermions transformed into a boson system by means of the boson 
expansion theory,^^ where the boson system is described by K-kinds of boson 
operators {Bj, Bj'; j=l,...K). 
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The main idea of Sakata^^ is that, just as in the classical theory a disso- 
lution of integrability (with the KAM mechanism) simply means the onset 
of chaotic motion/'^ in quantum systems a dissolution of quantum numbers 
may indicate the onset of quantum chaos. In accordance with the above 
definition of quantum chaos, we may classify the eigenstates \i > of the 
many-body Hamiltonian into three characteristic cases with the aid of the 
(/i, z/)-basis states, defined as |yU, z/ >= \fii...fiL,i^L+i---t^K >, which are spec- 
ified by /f-kinds of quantum numbers. 

1) Quantum integrable states : If one finds one of the |/i, u >-basis states for 
a given eigenstate \i > satisfying |</i, z/|i>p = l, then |i > is classified as 
a quantum integrable state, see Figure 8(a). 

2) Quantum KAM states : If \i > is described perturbatively starting from 
the l/i, z/ >-basis state, then it is a quantum KAM state, see Figure 8(b). 

3) Quantum chaotic states : If |i > is not described perturbatively starting 
from the |/i, u >-basis state, then it is regarded as a quantum chaotic state, 
see Figure 8(c). 

For further discussion of Sakata's approach see Ref. 35 and references 
quoted therein. 

Recently, another approach to the order-chaos transition has been pro- 
posed by Soloviev in the framework of nuclear structure. ^^'^^ Soloviev's main 
idea was to discuss the order to chaos transition in terms of the properties 
of nuclear wave- functions and to analyze how the structure of nuclear states 
changes with increasing excitation energy. He focused his attention on non- 
rotational states of rigid nuclei. The main conclusions are the following: 
1) Qrder is governed by the large components of the wave function of the 
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excited states. 

2) Chaos takes place in the small components of the wave function of the nu- 
clear excited states. The excited states are chaotic if their wave functions are 
composed of only small components of many-quasiparticle or many-phonon 
configurations. 

In our opinion these two approaches are quite similar and the Sakata 
approach also gives a simpler picture of the three cases discussed (regular, 
KAM, chaotic). 

6 Dynamical Quantum Chaos and Time Scales 

In this section we analyze in more detail the previously discussed problems 
of the chaotic time evolution for a quantum system. 

A quantum system evolves according to a linear equation and this is an 
important feature which makes it different from a classical system, for which 
the equations of motion can be nonlinear. On the other hand, the Liou- 
ville equation of the density function is linear both in classical and quantum 
mechanics but the evolution operator of Liouville has different spectral prop- 
erties. The classical Liouville operator has a continuous spectrum and this 
implies and allows chaotic motion. Instead, for bound systems, the quantum 
Liouville operator has a purely discrete spectrum, therefore no long-term 
chaotic behaviour. ^'^ 

As shown by Casati, Chirikov and coworkers^^'^° studying toy models like 
the kicked rotor, the time evolution of a quantum state follows the classical 
one, including the chaotic phenomena, up to a break time Ib- After that, in 
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contrast to classical dynamics, we get localization (dynamical localization). 
This means that persistent chaotic behaviour in the evolution of the quantum 
states and observables is not possible. Roughly speaking, chaotic behaviour 
is possible in quantum mechanics only as a transient with lifetime ts- The 
phenomenon of localization is clearly illustrated in Figure 9. This plot shows, 
in the case of the so-called standard map, the classical (solid curve) and 
quantum (dotted curve) unperturbed energy as a function of r (number of 
map iteractions). 

The value of the break time ts depends on the model studied and its exact 
behaviour is still controversial, but can be estimated from the Heisenberg 
indetermination principle as 

h 

where AE is the mean spacing of energy levels. The discrete spectrum of 
the Liouville operator cannot be resolved if t < t^, i.e. ts is the time at 
which the quantal evolution (of a wave packet, for example) "realizes" that 
the spectrum of the evolution operator is discrete. According to the Thomas- 
Fermi rule, AE ex h , where N is the number of degrees of freedom, i.e. the 
dimension of the configuration space. So, as h ^ 0, the break time diverges 
as tB ~ h^^ , and it does so faster, the higher A^ is. 

In the case of classical chaos, another time scale, the random time scale 
tji, much shorter than ts, can be introduced to estimate the time at which 
classical exponential spreading reaches the quantal resolution of the phase 
space. Thus, it is the full time for the exponential spreading of the minimum 
initial wave packet. As shown by Berman and Zaslavski^^ the random time 
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(called the breaking time by the two authors) follows the logarithmic law 

tn ^ A-i In (|) , (26) 

where S* is a classical action and A is the Lyapunov exponent of the sys- 
tem. In particular, Berman and Zaslavski have studied two different models, 
which allow one to calculate the random time scale t^. The first model is a 
periodically kicked oscillator, which gives t^j = ln{cost/h). The second is an 
ensemble of A^ atoms interacting with light in the resonant cavity. For this 
model tfl = In {cost N). 

6.1 Mean-Field Approximation and Dynamical Chaos 

Let us consider a A^-body quantum system with Hamiltonian H. The ex- 
act time- dependent Schrodinger equation can be obtained by imposing the 
quantum last action principle on the Dirac action 

S = j dt< ^{t)\m- - H\m > , (27) 

where ip is the many-body wavefunction of the system.^^ Looking for station- 
ary points of S with respect to variation of the conjugate wavefunction ip* 
gives 

th^^ij = Hij . (28) 

As is well known, it is usually impossible to obtain the exact solution of the 
many-body Schrodinger equation and some approximation must be used. 
In the mean-field approximation the total wavefunction is assumed to be 
composed of independent particles, i.e. it can be written as a product of 
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single-particle wavefunctions 0^. In the case of identical fermions, ip must be 
antisymmetrized. By looking for stationary action with respect to variation 
of a particular single-particle conjugate wavefunction 0* one finds a time- 
dependent Hartree-Fock equation for each 0^: 

where /i is a one-body operator. ^^ The main point is that, in general, the one- 
body operator h is nonlinear. Thus the Hartree-Fock equations are non-linear 
(integro-)differential equations. These equations can give rise, in some cases, 
to chaotic behaviour (dynamical chaos) of the mean-field wavefunction. 

In the mean-field approximation the mathematical origin of dynamical 
chaos resides in the nonlinearity of the Hartree-Fock equations. These equa- 
tions provide an approximate description, the best independent-particle de- 
scription, which describes, for a certain time interval, the very complicated 
evolution of the true many-body system. Two questions then arise: 

1) Does this chaotic behavior persist in time? 

2) What is the best physical situation to observe this kind of nonlinearity? 

To answer the first question, it should be stressed that, as shown pre- 
viously, quantum systems evolve according to a linear equation. Since the 
Schrodinger equation is linear, so is any of its projections. Its time evolution 
follows the classical one, including chaotic behaviour, up to the break time 
Ib- After that, in contrast to the classical dynamics, we get dynamical lo- 
calization. This means that persistent chaotic behaviour in the evolution of 
the states and observables is not possible. Nevertheless, we have seen that 
ts ^ h ~ , where A^ is the number of degrees of freedom of the system, thus 
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for a large number of particles the break time can be very long. 

Concerning the second question, it is useful to remember that, in the 
thermodynamic limit, i.e. when the number A^ of particles tends to infinity at 
constant density, the energy spectrum of the system is, in general, continuous 
and true chaotic phenomena are not excluded. ^^ 

When the mean-field theory is a good approximation of the exact many- 
body problem, one can use the nonlinear mean-field equations to estimate 
the transient chaotic behaviour of the many-body system. An important case 
where such an approach could be applied is the Bose-Einstein condensate of 
weakly-interacting alkali-metal atoms. ^^ In particular, we suggest that the 
collective oscillations of the Bose condensate^^'^^ can give rise to dynamical 
chaos. 

7 Conclusions 

In this paper we have reviewed various definitions of quantum chaos. In our 
opinion the simplest and clearest approach is that of Quantum Chaology, 
i.e. the study of quantum systems which are classically chaotic. Obviously, 
Quantum Chaology has some limitations, mainly it excludes systems without 
a classical analog, but its predictions in connection to the Random Matrix 
Theory are very accurate. 

Nowadays there are at least two important problems under investigation: 
the study of chaos without classical analog and the transient chaoticity of 
quantum systems. The problem of chaos in systems without a clear classical 
analog is very intricate and new ideas, like those of Sakata and Soloviev 
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discussed here, are needed. Also the precise behaviour of the time scales 
of dynamical chaos in quantum systems is not fully understood but some 
remarks can be made. 

We observe that the limitation to persistent chaotic dynamics in quantum 
systems does not apply if the spectrum of the Hamiltonian operator H is 
continuous. In the thermodynamic limit, i.e. when the number A^ of particles 
tends to infinity at constant density, the spectrum is, in general, continuous 
and true chaotic phenomena are not excluded. We have seen that the break 
time ts is very long for systems with many particles. The transient chaotic 
dynamics of quantum states and observables can be experimentally observed 
in many-body quantum systems. Moreover, the fact that the break time ts 
increases with the number of microscopic degrees of freedom explains the 
chaotic behaviour of macroscopic systems, without invoking a role for the 
observer or the environment. 

The study of quantum chaos for many-body or continuum systems (field 
theory) is a very promising field of research which can also help to better 
understand the foundations of quantum theory and statistical mechanics. 

* * * 

One of us (L.S.) is greatly indebted to Prof. M. Robnik for many sugges- 
tions. 
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Figure Captions 

Figure 1: P{s) for "cold" deformed rare-earth nuclei (adapted from Ref. 

17). 

Figure 2: Spectral statistics for nuclei with atomic mass 24 < A < 244 

and excitation energy of few MeV. a) 2+ and 4+ states, even-even nuclei; b) 

all other states, even-even nuclei; c) states with non-natural parity, odd-odd 

nuclei; d) states with natural parity, odd nuclei (adapted from Ref. 17). 

Figure 3: Comparison of the nearest-neighbour spacing distribution P{s) 

and spectral rigidity A3(L) of the Nuclear Data Ensemble (NDE) with the 

GOE predictions (adapted from Ref. 30). 

Figure 4: Histograms of the level distances of quantum energies in the 

Hydrogen atom in magnetic fields at different values of the scaled energy. 

The smooth curves are the results of the fits to the histograms (adapted 

from Ref. 32). 

Figure 5: Spectral rigidity A3(L) for energy level sequences of the Hydrogen 

atom in a magnetic field of 6 Tesla in three different energy intervals. The 

transition to the GOE distribution as soon as the classical motion becomes 

chaotic is also visible here (adapted from Ref. 32). 

Figure 6: The Poincare sections of the model. From the top: f = 1, f = 1.1 

and V = 1.2. Energy E = 10 and interaction g = 1 (adapted from Ref. 28). 

Figure 7: P{s) distribution. From the top: v = 1 {u = 0.92), v = 1.1 

{uj = 0.34) and v = 1.2 {u = 0.01), where u is the Brody parameter. First 

100 energy levels and interaction g = 1. The dotted, dashed and solid curves 

stand for Wigner, Poisson and Brody distributions, respectively (adapted 
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from Ref. 28). 

Figure 8: Exact eigenstates \i > expressed in tlie |/i, z/ > basis states: a) 

Quantum Integrable States; b) Quantum KAM States; c) Quantum Cliaotic 

States (adapted from Ref. 35). 

Figure 9: Classical (solid curve) and quantum (dotted curve) unperturbed 

energy < n^(r) >= 2E as a function of time r (number of map iterations) 

(adapted from Ref. 5). 
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